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Abstract 

The six-fermion production processes e'’"e“ —> with all 

the lepton flavours and q = u,d, c, s, relevant to the study of the 
intermediate-mass Higgs boson at future e^e“ linear colliders, are 
analysed. A Monte Carlo program, taking into account the whole set 
of tree-level scattering amplitudes and the relevant radiative effects, 
is developed to provide integrated cross sections and generation of 
unweighted events. The complete calculation is compared with the 
available results of real Higgs production, and the opportunities of 
precision studies with event generation are discussed, demonstrating 
the relevance of a full six-fermion calculation. Numerical results for 
integrated cross sections with various kinematical cuts and including 
radiative effects are given and commented. In the analysis of event 
samples, several distributions are studied and found to be sensitive to 
the presence and to the properties of the Higgs boson. 
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1 Introduction 


The search for the Higgs boson and the stndy of its properties will be among 
the most important tasks of elementary particle physics at future e’''e“ linear 
colliders at the TeV scale (NLC) [^. 

The research carried on at present colliders can explore an interval of 
Higgs masses below ~ 100 GeV at LEP2 0 or possibly 120-130 GeV at the 
upgraded Tevatron |0]. The remaining mass range, up to ~ 800 GeV, will 
be in the reach of the future colliders LHG Q and NLG. In particular the 
precision studies that will be possible in the clean environment of NLG will 
be of great help in the determination of the Higgs boson properties. 

A range of particular interest for the Higgs mass is between 100 and 200 
GeV, as many arguments both of experimental and of theoretical nature 
indicate. Indeed a lower limit of ~ 90 GeV is given by recent results in the 
direct search at LEP2 p|, while from fits to electroweak precision data an 
upper limit of ~ 280 GeV at 95% G.L. is obtained p. 

In this range two mass intervals may be considered: for mn < 140 GeV 
the Higgs decays mainly into bb pairs, while for thh > 140 GeV the decays 
into WW and ZZ pairs become dominant. Therefore in the first case the 
mechanisms of Higgs production relevant to e+e“ colliders, Higgs-strahlung 
and VViy = IT, Z) fusion, give rise to signatures that contain four fermions 
in the final state, which have been extensively studied in the recent past |^- 
[^]. In the second case, in which mn > 140 GeV, six fermions are produced in 
the final state. More generally, six-fermion (6/) final states come from other 
relevant processes at NLG, such as ti and three-gauge-boson production [PI]|] . 

As shown by complete four-fermion calculations for WW and light Higgs 
physics at LEP2 i-i 0, 0, full calculations of 6/ production processes 
allow one to keep phenomenologically relevant issues under control, such 
as off-shellness and interference effects, background contributions and spin 
correlations. Some calculations of such processes e+e“ —> 6f have recently 
been performed |]0]-|0], with regard to fop-quark, Higgs boson and WWZ 


physics at NLG. Moreover recent progress in the calculation of processes 
> 6 jets [0 and of 2 —> up to 8 partons QGD amplitudes ||0, 0 


e^e 


should be mentioned for their relevance in QGD tests at lepton and hadron 
machines. These calculations rely upon different computational techniques, 
such as helicity amplitude methods for the evaluation of the very large num¬ 
ber of Feynman diagrams associated to the process under examination, or 
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iterative numerical algorithms, where the transition amplitudes are calcu¬ 
lated numerically without using Feynman diagrams. 

Concerning Higgs physics, an analysis of the processes e+e“ —>• fi~^n~udT~9T- 
and e+e“ — ^~ude~ Ve has been performed in ref. [|T^, where the Higgs can 
be produced by Higgs-strahlung and the subsequent decay proceeds through 
W^W~ pairs. Special attention has been devoted to the calculation of the 
Higgs boson signal and of its Standard Model background, with special em¬ 
phasis on the determination and analysis of angular correlation variables, 
particularly sensitive to the presence and to the spinless nature of the Higgs 
particle. The 6/ hnal states, where the Higgs signal gives rise to two charged 
currents, have also been considered in ref. |^, studying cross sections and 
invariant mass distributions for the processes e+e“ —>• ffqq'f f"- 

The case of the Higgs boson decay in neutral currents has been briefly 
addressed for the signal alone in ref. [|^ with the study of the reaction 
—>• e'^e~i'e,VeUU. The aim of the present paper is to complete and ex¬ 
tend the analysis of ref. to include general qq neutral currents contribu¬ 
tions and the effect of the contributions from undetectable different-flavour 
neutrinos, in such a way as to provide realistic predictions for processes 
e’''e“ —>■ l^l~h>9qq at the parton level. In the following, 6-quark tagging will 
be assumed, leaving aside bb hnal states, which lead to an interplay between 
Higgs and top physics and will be studied elsewhere. Consisting of only 
two jets, the processes considered in the present paper are free from QCD 
backgrounds. 

The outline of the paper is as follows. In Section 2 the physical process 
is presented and the main technical issues of the calculation are explained. 

In Section 3 several numerical results are shown and discussed and the po¬ 
tentials of full 6/ calculations are stressed. Finally, Section 4 contains the 
conclusions. 


2 Physical process and computing technique 

The production of an intermediate mass Higgs boson is studied in the process 
e’''e“ —>■ qql~^l~uV, where a sum is made over the contributions of the u, d, c 
and s quarks, of the three neutrinos and of / = e, /r, r. Particular attention 
will be devoted to the signature qqe^e~i'V. 

One of the interesting features of this process is the presence of both 
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charged current and neutral current contributions |]T3, which is a situation 


never studied before, since all the six-fermion signatures analysed in the lit¬ 
erature [0-[0 involve only charged current contributions. Moreover this 


class of processes receives contribution from diagrams with up to three t- 
channel gauge bosons. This feature is of particular interest because of the 
large centre-of-mass (c.m.) energy, -^/s, at which the NLC will operate. These 
topologies enhance the cross section with growing s. The capability to pro¬ 
vide predictions for processes with many electrons and electron-neutrinos 
in the final state is therefore crucial to discuss NLC physics. The present 
study demonstrates the possibility to deal successfully with the dynamics 
calculation and phase-space integration of this class of final states. Another 
important property is that the process is free from QCD backgrounds be¬ 
cause only two jets are produced. As a drawback, the total cross section is 
smaller than in the 6/ processes with four or six jets. 

However, the sums over quark, charged lepton and neutrino flavours, as 
well as the combined action of different mechanisms of production (see £g. |l|), 
contribute to give a signihcant size to the cross section, so that, assuming 
an integrated luminosity of 500 fb“^/yr and a Higgs mass of, say, 185 GeV, 
more than 1000 events can be expected at a c.m. energy of 360 GeV and 
more than 2000 at 800 GeV (see fig. |^). In particular, as will be seen in the 
numerical results, the presence of the t-channel contributions of vector boson 
fusion gives an enhancement of the cross sections at very high energies. 

The diagrams containing a resonant Higgs boson coupled to gauge bosons 
for the qqe^e~vV hnal state are shown in fig. |I[ As can be seen, there are four 
terms of the Higgs-strahlung type and two of the fusion type. At relatively 
low energies, ^/s < 500 GeV, the process of Higgs-strahlung dominates and, 
in particular, the charged current term is the most important one. As the 
energy is increased, the t-channel terms of vector boson fusion become more 
and more important, as they grow with increasing s, and they are dominant 
above 500 GeV. The diagrams for the processes with or instead 
of e’''e“ and/or instead of z/et'e in the hnal state are a subset of those 

illustrated here. 


The full set of diagrams containing a physical Higgs boson coupled to 
gauge bosons includes also those shown in £g. ||. These contributions are 
non-resonant, as the Higgs is exchanged in the t-channel, and their size can 
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Figure 1: Feynman diagrams for the process e'’"e —> qqe'^e vv with a resonant 

Higgs boson. 

be expected to be negligible at low energies; at high energies, however, they 
play an important role in preserving gauge invariance and unitarity of the 
S'-matrix. This point will be discussed in more detail in the next section. 


The total number of tree-level Feynman diagrams for the process under 
examination is of the order of one thousand, which makes the calculations 
very complicated for what concerns the determination of the transition ma¬ 
trix element as well as the phase-space integration. 

The transition matrix element is calculated by means of ALPHA [p0|] , 
an algorithm that allows the calculation of scattering amplitudes at the tree 
level without the use of Feynman diagrams. The results produced by this 
algorithm in a large number of calculations of multi-particle production are in 
full agreement with those obtained by programs using traditional techniques 
This fact may be considered as a signihcant test of ALPHA. 


12, 20-23 


Anyway, also in the present work, checks have been made, reproducing 
by means of the helicity amplitude method p3 some of the results given by 
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Figure 2: Feynman diagrams with a non-resonant Higgs boson (two similar dia¬ 
grams are also obtained from these by particle exchanges). 


ALPHA for the Higgs “signal” (which is dehned below) and hnding perfect 
agreement. 

For the integration over the phase space, as was already done in refs. 
[pAl , [T^, a code has been developed by adapting the Monte Carlo program 
HIGGSPV/WWGENPV [^, originally developed to treat four-fermion 
production, to make 6/ calculations. The code can be used to perform Monte 
Garlo integrations and obtain cross sections, or to generate samples of un¬ 
weighted events. Kinematical cuts can be introduced to simulate realistic 
experimental conditions. The effects of initial-state radiation (ISR) and 
beamstrahlung are taken into account by means of the standard convo¬ 
lution formula 


cr = J dzidz 2 DBs{zi,Z 2 ;s) J dxidx 2 D{xi,s)D{x 2 ,s)d{zi,Z 2 ]Xi,X 2 ;s) . (1) 

An accurate importance sampling procedure is required in the Monte 
Garlo integration to take care of the complicated structure of “singularities” 
in the integrand. This structure results from the presence of several mech¬ 
anisms of Higgs production, and also from additional sources of variance 
among the very large number of background diagrams present in the matrix 
element. 

The “singularities” given by different terms correspond to different regions 
of the (14-dimensional) phase space and in general must be treated with 
different sets of integration variables. As a consequence, a multichannel 
importance sampling technique is needed. If n channels are introduced, the 
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Figure 3: s-channel and t-channel topologies considered in the importance sam¬ 
pling. 

integral is written as 

7X ^ \ Tl 

j /(x)rf/i(x) = ^j p(x) = ^Pi(x), (2) 

where each x^*) is a set of integration variables with a corresponding measure 
dfii, and pi is a suitably normalized distribution function for the importance 
sampling in the Fth channel. 

The choice of integration variables is made within two kinds of phase- 
space decompositions, corresponding to two diagram topologies: s-channel, 
based on the Higgs-strahlung terms, and f-channel, based on the fusion terms, 
as illustrated in £g. 

For the s-channel topology the phase-space decomposition reads: 

d^Q{P; qi, ..., qq) = (27r)^^d<F2(-P; Qi, Q2)d^2iQi] Qi, Q2) 

d^2{Q2] Q 3 , Qi)d^2{Q3] ^35 (li)d^2{Q4] ? 5 ) Qg) 
dQldQldQldQl (3) 

where P is the initial total momentum, qi are the momenta of the out¬ 
going particles, while Qi are those of the internal particles. The notation 
d^2{QuQj,Qk) indicates the two-particle phase space; the momenta of the 
hnal particles Qj, Qk are first generated in the rest frame of the particle of 
momentum Qi and then boosted to the laboratory frame. The integration 
variables are the four invariant masses Q^,... ,Q\ and five pairs of angular 
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variables cos0j, (pi, one for each <^$2 term. The invariant masses are sampled 
according to Breit-Wigner distributions of the form 


(M2 - Q2)2 + V y 

given by the propagators of the Higgs or gauge bosons in the internal lines 
{N is a normalization factor). For the angular variables a flat distribution 
is assumed. The various s-channel terms differ for permutations of the ex¬ 
ternal momenta and for the parameters T, M in the importance sampling 
distributions. 

In the case of f-channel diagrams, the phase space is 

d^e{P-, gi,..., qe) = {2Tr)^d^3{P; Qi, qi, ga) 

d^ 2 {Qi] Q 2 , Q3)d^2{Q2] gsj qi)d^2{Q3] q^-, q%) 
dQldQldQl ( 5 ) 


where, as before, the g* are the outgoing momenta, while the Qi are internal 
time-like momenta. The integration variables are the three invariant masses, 
Q\, Q 2 , Qi, three pairs of angular variables cos 6 *, 0 relative to the three d ^2 
terms, and, for the three-body phase space d^s, one energy, g°, and four 
angular variables, cos 6 * 1 , 0i, cos 6 * 2 , 02 - The invariant masses are sampled, as 
in the s-channel case, according to Breit-Wigner distributions; one angular 
variable in ^$ 3 , say cos^i, is sampled by means of the distribution 


(M2 +^gO(l-COS0i))2 + r2.M2’ ' ' 

corresponding to the propagator of one space-like gauge boson (V = W, Z) 
emitted by the initial electron or positron (typically one of the bosons par¬ 
ticipating in the fusion into Higgs); in some channels, corresponding to back¬ 
ground diagrams, also another angular variable cos 0 relative to a two-body 
term d^ 2 , is sampled in a similar way, in order to take into account the “sin¬ 
gularity” associated with a boson propagator. All other variables have flat 
distributions. 


3 Numerical results and discussion 


The numerical results presented in this section are obtained with the same set 
of phenomenological parameters as adopted in ref. [O. Namely, the input 
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parameters are Mw and and other qnantities, snch as sin^ 9w, « 
and the widths of the W and Z bosons, are compnted at tree level in terms 
of these constants. The Higgs width inclndes the fermionic contribntions 
h —> /ip, rr, cc, bb, with running masses for the quarks (to take into account 
QCD corrections [§), the gluonic contribution h ^ gg [§, and the two- 
vector boson channel, according to ref. [^. The denominators of the bosonic 


propagators are of the form +iTM, with hxed widths T. As already 

discussed in ref. [|T^, the aim of this choice is to minimize the possible sources 
of gauge violation in the computation . 

Such gauge violations have been studied by the same methods as were 
used in ref. |]^. In particular, for what concerns SU{2) gauge symmetry. 


comparisons have been made with results in the so-called “fudge scheme” [pi|] . 
A disagreement has been found at the level of few per cent for a c.m. energy of 
360 GeV. The disagreement vanishes at higher energies. By careful inspection 
of the various contributions, it has been checked that the deviation at lower 
energies is due to the well-known fact that a given fudge factor, close to 
a resonance, mistreats the contributions that do not resonate in the same 
channel. Concerning U{1) invariance, a test has been performed by using 
different forms of the photon propagator and hnding perfect agreement, up 
to numerical precision, in the values of the squared matrix element. 

The hrst group of results discussed in this section refers to cross-section 
calculations, performed by using the program as an integrator of weighted 
events. The signature considered in the hrst plots of total cross section is 
qql~^l~uV, where, in addition to the sums over quark and neutrinos havours 
already mentioned, there is a sum over I = e, p, r. All other results are 
instead relative to the signature qqe^e~vV. Some samples of unweighted 
events, obtained by using the code as a generator, are then analysed in the 
remaining part of this section. 


3.1 Total cross sections 

In £g. ^ the total cross section (including the contribution of all the tree- 
level Feynman diagrams) is shown for three values of the Higgs mass in the 
intermediate range, 165,185 and 255 GeV, at energies between 360 and 800 
GeV. To make a hrst analysis, the following kinematical cuts are adopted: 
the invariant mass of the quark-antiquark pair and that of the charged lepton 
pair are required to be greater than 70 GeV, the angles of the two charged 
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Figure 4: Total cross section for the process e'’'e“ —> qqyi~uV in the Born 
approximation, as a function of ^/s for three different values of the Higgs mass 
iriH- The angles 9{l~) of the charged leptons with the beam axis are in the 

interval 5°-175°, the e^e“ and the qq invariant masses are larger than 70 GeV. 


leptons with respect to the beam axis within 5° and 175°. This choice is 
applied to the quantities shown in hgs. and 1^. Another set of cuts, with 
a lower limit of 20 GeV on the invariant mass, is adopted in £g. ^ and 
in the study of event samples. 

The increase with energy, common to all three curves in £g. |^, is due, at 
high energies, to the f-channel contributions; in the case of rriH = 255 GeV, 
the steep rise near ^/s = 360 GeV is related to the existence of a threshold 
effect for the Higgs-strahlung process at an energy y/s ~ tuh + Mz- 


In fig. ^ the total cross section is plotted, with the cut on the invariant 
mass of the charged lepton pair reduced to 20 GeV. The effect of this modifi¬ 
cation, as expected, is an enhancement of the cross section in the low-energy 
region: indeed the most important contribution at energies below 500 GeV is 
given by the Higgs-strahlung diagram, with the Higgs decaying into two W 
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Figure 5: The same as in fig. ^ with the cut on the l^l invariant mass reduced 
to 20 GeV. 

bosons, which will be indicated from now on as the charged-current Higgs- 
strahlung diagram, and which is characterized by a broad distribution of the 
invariant mass that goes well below 70 GeV. This cut is still sufficient 
to reduce to a negligible size the contribution of virtual photon conversion 
into l^l~ pairs. The behaviour of the cross section as the Higgs mass is var¬ 
ied depends on the interplay of the various production mechanisms and of 
the decay branchings involved; this behaviour can be better observed in the 
“signal” contribution that will be dehned below (see £g. 

The effect of ISR is illustrated in £g. ||, for a Higgs mass of 185 GeV and 
for the signature qqe~^e~uV, to which all the remaining results refer. Here 
the e’''e“ invariant mass is again greater than 70 GeV. 

The cross section in the presence of ISR is lowered by a quantity of the 
order of 10% with respect to the Born approximation. This phenomenon 
can be easily understood, since the initial-state radiation reduces the c.m. 
energy, so in this case it produces a shift towards energy values where the 
cross section is smaller. 
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Figure 6: Effect of initial-state radiation on the total cross section of the process 
e+e“ ^ qqe^e~vV as a function of \fs for a Higgs mass mu = 185 GeV. Cuts are 
the same as in fig. 

3.2 Definition and study of the Higgs signal 

The results discussed so far, as stated above, are given by the sum of all the 
tree-level Feynman diagrams. Strictly speaking, this is the only meaningful 
procedure. On the other hand, there is a number of reasons to consider a 
subset of diagrams that can be dehned as the Higgs signal and to dehne a 
corresponding background. In the hrst place this is of great interest from the 
point of view of the search for the Higgs boson in the experiments. More¬ 
over, as will be shown, such a dehnition allows one to make a comparison 
with results obtained in the narrow width approximation (NWA) [0^ - 
which are the only available estimations unless a complete 6/ calculation is 
performed. In principle, whenever a subset of diagrams is singled out, gauge 
invariance may be lost and unitarity problems may arise. However, in the 
following, an operative dehnition of signal and background is considered and 
its reliability is studied for various Higgs masses and c.m. energies. 

The signal is dehned as the sum of the six graphs containing a resonant 
Higgs boson, shown in hg. |l|. The background is dehned as the sum of all 
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the diagrams without a Higgs boson. In this dehnition the diagrams with 
a non-resonant Higgs boson coupled to gauge bosons, shown in £g. ||, are 
missing both in the signal and in the background. Such a choice has been 
dictated by the fact that these non-resonant contributions cannot correctly 
be included in the signal, since they cannot find a counterpart in the NWA, 
and because of gauge cancellations with background contributions at high 
energies; however, as they depend on the Higgs mass, they should not be 
included in the background as well. 

In order to give a quantitative estimate of the validity of this dehnition, 
the total cross section (sum of all the tree-level 6/ Feynman diagrams) is 
compared in £g. with the incoherent sum of signal and background. The 
cuts are as in kg. in particular with the e’''e“ invariant mass greater than 
70 GeV. In order to understand the meaning of these results, it is impor¬ 
tant to note that, as observed above, the contributions of the diagrams of 
£g. 1^ are absent both from the signal and from the background: thus if we 
indicate these contributions to the scattering amplitude as A^t, and the sig¬ 
nal and background amplitudes as Ag and A^ respectively, the total squared 
amplitude is 

{A]"^ = \As + Ab + Aht\‘^ . (7) 

The terms neglected in the incoherent sum of signal and background are | Aht \ ^ 
and all the interference terms. Among these, the interferences of Aht with 
the rest are dominant at high energies as they involve gauge cancellations. 


The curves of £g. ^ show that up to 500 GeV the total cross section and 
the incoherent sum of signal and background are indistinguishable at a level 
of accuracy of 1%, and the dehnition of signal may be considered meaningful; 
at higher energies, this separation of signal and background starts to be less 
reliable, since it requires to neglect ehects that are relevant at this accuracy. 
In particular, at 800 GeV the deviation is of the order of a few per cent and 
it decreases when the Higgs mass passes from 165 to 185 and to 255 GeV. 

The above results are obtained with the set of kinematical cuts in which 
the e’''e“ invariant mass is greater than 70 GeV, but when this cut is at 20 
GeV, the diherence between full cross section and incoherent sum of signal 
and background is signihcantly reduced (about 3-4% at 800 GeV). The anal¬ 
ysis of event samples presented in the following is made within this latter set 
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Figure 7: Full six-fermion cross section compared with the incoherent sum of “sig¬ 
nal” (S) and “background” (B). A detailed discussion and an operative definition 
of “signal” and “background” are given in the text. 
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of cuts, so that, up to 800 GeV, it can be considered reliable, at the level of 
accuracy of a few per cent, to speak of “background”, as will be done. 

On the other hand the problems arising when a dehnition of signal and 
background is attempted show the importance of a calculation involving the 
full set of tree-level Feynman diagrams to obtain reliable results, especially 
at high energies. 

A comparison with the NWA is shown in £g. ^ for the processes e’''e“ —> 
qqe^e~vV and e’''e“ —*■ qq^~^fi~i>V, where, for the sake of comparison, no 
kinematical cuts are applied and the results are in the Born approximation. 
Here agig is the signal cross section, containing the contributions of the six 
diagrams of £g. |I] (or the suitable subset of these for the case of the hnal state 
qqfi ~^and their interferences. The cross section in the NWA, is 

obtained in the following way (for dehniteness the case with e’''e“ in the hnal 
state is considered): the known cross sections for the processes of real Higgs 
production —>• huV, he^e~ [^, and Zh [^] are multiplied 


by the appropriate branching ratios, so as to obtain six terms corresponding 
to the diagrams of kg. |^; then the incoherent sum of these terms is taken. 
Thus the comparison between agig and aNWA gives a measure of interfer¬ 
ence between the different production mechanisms and of off-shellness effects 
together. 


As can be seen in kg. the relative dikerence R is of the order of some 
per cent, depending on the Higgs mass and the c.m. energy; in some cases 
it reaches values of more than 10%, with no substantial difference between 
the two final states considered. In particular the ok-shellness effects are 
much more important than the interference ones. In fact the relative size 
of the interferences has been separately evaluated by means of a comparison 
between agig and the incoherent sum of the six diagrams of kg. ^ and has 
been found to be at most 2%, but generally less than 1% for the c.m. energies 
and Higgs masses considered here. 

The size of the ok-shellness ekects, comparable with the ISR lowering, 
indicates the importance of a full 6/ calculation in order to obtain sensible 
phenomenological predictions. 

In kg. ^ the signal cross section is shown as a function of the Higgs mass 
for dikerent c.m. energies. The behaviour is related to the branching ratios of 
the decays of the Higgs boson into gauge bosons and the dikerences between 
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Figure 8: Comparison between the signal cross section obtained by a diagram¬ 
matic six-fermion calculation and the one calculated in the narrow width approx¬ 
imation (see the discussion in the text), as a function of ^/s (upper row) and of 
the Higgs mass (lower row). 
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Figure 9: Signal cross section as a function of the Higgs mass for three different 
c.m. energies. 


the three energy values considered are due to the variations in the relative 
sizes of the different signal contributions at different energies. 


3.3 Distributions 

The results presented in the following refer to samples of unweighted events 
for a Higgs mass of 185 GeV at energies of 360 and 800 GeV, with or without 
the effect of ISR and beamstrahlung. The cuts adopted in all cases are the 
following: the invariant mass of the qq pair greater than 70 GeV, the invariant 
mass of the e’''e“ pair greater than 20 GeV, and the angles of the electron 
and positron with respect to the beam axis between 5° and 175°; further 
cuts, applied in the analysis of particular cases, will be described later. The 
numbers of events in all the samples are normalized to the same luminosity. 

In £g. |T^ the invariant masses of two different systems of four momenta 
are studied at c.m. energies of 360 and 800 GeV in the Born approxima¬ 
tion (dashed histograms) as well as with ISR and beamstrahlung (solid his- 
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tograms). The first set is given by e+e“+ missing momentum, where the 
missing momentum is dehned as qmiss = pfn + Pin ~ Qe+ ~ <?e- ~ Qq ~ Qq- 
In the Born approximation this set of momenta corresponds to the system 
e^e~h'V. The other set considered is that corresponding to the four-fermion 
system qqe^e~. As can be seen in £g. |I|, these are two of the possible sets 
of four fermions produced by the decay of the Higgs boson in the process 
under consideration; there is also a third set, qqvV, whose invariant mass 
distribution, however, does not contain any new feature. The presence of the 
Higgs boson can be revealed by a peak in the distributions of these invariant 
masses. Indeed, in the Born approximation (dashed histograms), a sharp 
peak around 185 GeV can be seen in each of the histograms of £g. At 
a c.m. energy of 360 GeV, the most remarkable one is that of e+e“-|- miss¬ 
ing momentum, associated to the system e+e“ w, as it receives contributions 
from the charged current Higgs-strahlung diagram, which is dominant at this 
energy. In the presence of ISR and beamstrahlung, this peak is considerably 
lowered and broadened, while the other distribution, not involving the miss¬ 
ing momentum, is not significantly affected by radiative effects. At 800 GeV 
this phenomenon is even more evident, because the peak in the first distri¬ 
bution is completely eliminated by radiative effects, as a consequence of the 
small size of the charged current Higgs-strahlung term at this energy, while 
the second distribution results to be very sensitive to the presence of the 
Higgs, since it receives, around 185 GeV, contributions from the diagram of 
WW fusion into Higgs, which is the dominant signal term at high energies, 
and the presence of ISR and beamstrahlung does not modify the shape of 
the histogram. Thus, at high energies, a very clean signal of the Higgs boson 
is provided by the process under study through this distribution. 


The quantities analysed above are useful to reveal the presence of the 
Higgs boson and to determine its mass. Other variables can be considered to 
study the properties of this particle, such as spin and parity. Some examples 
are considered in figures 00 and |T^. 

When the process e’''e“ ^ HZ is considered, a variable that can give 
evidence of the scalar nature of the Higgs is the angle 6z of the Z parti¬ 
cle direction with respect to the beam in the laboratory frame. It is well 
known that the differential cross section da/d cos 6z goes as sin^^^ at 


energies much greater than Mz and away from the threshold for Higgs pro¬ 
duction. A similar situation is expected to occur for the 6/ process under 
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rrih = 185 GcV (Born vs. ISR + BS) 






Figure 10; Invariant-mass distributions for four-fermion systems in the Born 
approximation (dashed histograms) and with ISR and beamstrahlung (solid his¬ 
tograms) at y/s = 360 GeV (upper row) and y/s = 800 GeV (lower row). 
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Figure 11: Distribution of the angle 9z of the qq pair with respect to the beam 
in the laboratory frame at ^/s = 360,800 GeV. The solid histogram represents the 
full calculation, the dashed histogram is the contribution of the background. 


study when the Higgs-strahlung contributions are dominant. The distribu¬ 
tion da/dOz is shown, at the c.m. energies of 360 and 800 GeV, in £g. 
where the Z particle is reconstructed as the sum of the quark and antiquark 
momenta (indeed this is the case for the dominant diagram). The contri¬ 
bution from the background alone (dashed histogram) is also shown. The 
shape of the solid histogram shows the expected behaviour at 360 GeV, where 
Higgs-strahlung dominates. 


At the c.m. energy of 800 GeV, where the dominant signal diagram 
is WW fusion into Higgs, the situation is substantially different, since the 
process of Higgs production is of the t-channel type. One variable that results 
very sensitive to the presence of the Higgs boson is shown in £g. 12 and 
indicated as cos 9zz', 9zz is the angle between the three-momenta in the 
laboratory frame of the qq and e’''e“ pairs, which correspond, in the diagram 
of WW fusion, to the Z particles coming from the Higgs. The full distribution 
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Figure 12: Distribution of the angle Ozz between the qq and the e^e“ pairs in 
the laboratory frame at \fs = 800 GeV. The solid histogram is the full calculation, 
the dashed histogram is the background. 

(solid line) and the contribution from the background alone (dashed line) are 
particularly distinguished in the region near 1. There is here a clear signal 
of the presence of the Higgs, and such a variable can be used to impose 
kinematical cuts to single out signal contributions. This phenomenon is 
however of a kinematical nature, and is not directly related to the scalar 
nature of the Higgs, but is rather a consequence of the smallness of diagrams 
with the same topology of the WW fusion in the background. 


Another variable has been considered at the energy of 800 GeV in £g. B 
It is the angle of the Z particle, reconstructed as the sum of the electron 
and positron momenta, with respect to the beam, in the rest frame of the 
system qqe^e~. The reference diagram is always the WW fusion: this may be 
regarded asymptotically as an s-channel WW scattering into ZZ, and, in the 
rest frame of the incoming WW pair, the angular distribution of the produced 
ZZ pair is determined by the scalar nature of the exchanged particle. In the 
hrst row of fig. the plot on the left is made without additional cuts. 
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while the plot on the right is obtained with the requirement cosOzz > 0, 
so as to reduce the background. In the second row the invariant mass of 
qqe^e~ is required to be smaller than 250 GeV (left) and within 20 GeV 
around the Higgs mass (right) in order to further suppress the background. 
A clear difference between the shape of the full distribution and that of the 
background can in fact be seen, and in the last three plots the behaviour 
is very similar to the sin^^ distribution expected on the basis of the above 
observations. 


4 Conclusions 

The processes e+e“ —>■ qql~^l~uT7 have been studied in connection with the 
search for an intermediate-mass Higgs boson. The study, which extends a 
previous analysis of 6/ signatures with only two jets, is characterized by the 
presence of neutral current contributions that were never considered before 
and by the fact that several mechanisms of Higgs production are simultane¬ 
ously active. 

The tool used for the numerical calculations is a Fortran code based on the 
algorithm ALPHA, for the determination of the scattering amplitude, and 
on a development of the Monte Garlo program HIGGSPV/WWGENPV, for 
the phase-space integration. 

The total cross section, including all the tree-level Feynman diagrams, 
has been calculated with various kinematical cuts and taking into account 
the effects of ISR and beamstrahlung. 

A dehnition of signal and background has been considered and its reli¬ 
ability has been studied. To this end the incoherent sum of “signal” and 
“background” has been compared with the full cross section, and this has 
shown deviations that, up to a c.m. energy of 500 GeV, are negligible to an 
accuracy of 1%, but may be of several per cent at 800 GeV (hg. These 
deviations are, however, reduced when the kinematical selection criteria be¬ 
come more inclusive. 

A comparison of the “signal” cross section with results in the NWA has 
shown that off-shellness effects have a relative size of several per cent (hg. |^). 

The results of hgs. ^ and || show the importance of a complete 6/ calcu¬ 
lation to produce reliable results at the TeV scale. 
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185 GeV Vs = 800 GeV 


m, = 




55 55 


Figure 13: Distribution of the angle of the e'’'e“ pair with respect to the 
beam axis in the rest frame of the qqe~^e~ system at ^/s = 800 GeV. The solid 
histogram is the full calculation, the dashed histogram is the background. First 
row: no additional cuts (left), cosOzz > 0 (right), where the angle 9zz is defined 
in the text and shown in fig. |T^. Second row: qqe'^e~ invariant mass smaller than 
250 GeV (left) and within 20 GeV around the Higgs mass (right). 
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In the study of generated events the problem of hnding observables that 
are sensitive to the presence of the Higgs and to its properties has been 
addressed. The presence of several mechanisms of Higgs production, whose 
relative importance varies with energy, requires that different variables be 
considered according to the energy range studied. 

The invariant masses of two sets of four fermions have been analysed hrst 
(fig- HI: one, relative to the system e^e +missing momentum, is relevant 
to the detection of the Higgs boson at 360 GeV of c.m. energy, but, at 
800 GeV, the effects of ISR and beamstrahlung prevent to study the Higgs 
by means of this distribution. The other invariant mass, relative to the 
system e^e~qq, is instead particularly useful at high energies and is almost 
completely unaffected by radiative effects. 

Three angular variables have then been studied: the angle 6z (fig. 0) 
is suited to reveal the spin zero nature of the Higgs at 360 GeV, where the 
Higgs-strahlung dominates, but it gives no information at 800 GeV. The 
angles 6zz (fig- [0) and 6*^ (£g. |T^ are very useful at 800 GeV: the first one 
is very effective to single out the signal, but is not able to distinguish the spin 
nature of the Higgs; the second one has a distribution whose shape is very 
different from that of the background, and is related to the spinless nature 
of the Higgs particle. 

The computing strategy and the relative computer code developed in this 
work have been applied to study intermediate-mass Higgs physics. However, 
the variety of diagram topologies present in the matrix element and taken into 
account in the Monte Garlo integration, as well as the possibility provided 
by ALPHA of dealing with any kind of process, including now also QGD 
amplitudes, give the opportunity to examine other topics relevant to physics 
at future colliders, where 6/ production is involved. 
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